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Abstract 



We consider an initially bound quantum particle subject to an external 

time-dependent field. When the external field is large, the particle shows a 

tendency to repeatedly return to its initial state, irrespective of whether the 

frequency of the field is sufficient for escape from the well. These recurrences, 

which are absent in a classical calculation, arise from the system evolving 

primarily like a free particle in the external field. 
PACS numbers: 42.50.Hz, 32.80.Fim 
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A very common experimental probe of a physical system is the study of its response 
to an imposed time-dependent perturbation. The behavior of quantum systems subject 
to strong time-dependent forces is of fundamental interest to a variety of situations such 
as field-enhanced decay in nuclei ]I], highly excited atoms in microwave fields || and rf- 
driven Josephson junctions ||. For weak external forces, a very successful approach is 
standard perturbation theory, where the response of a physical system is calculated in terms 
of correlation functions evaluated in the absence of the external perturbation. However, 
there are a steadily increasing number of experimental situations where the probe is so 
invasive that one does not expect perturbative approaches to provide even a qualitative 
picture of the behavior of the system. For example, the use of high intensity lasers in 
a variety of experiments during the last few years has revealed novel phenomena like the 
decreasing of ionization probability with increasing laser intensity |4j]. In these situations, 
there is no simple theoretical picture akin to perturbation theory and so the analysis of 
quantum systems under strong driving remains an unsolved problem. 

We consider the following problem: a quantum particle is in the ground state of a 
Hamiltonian whose spectrum is partly discrete and partly continuous. It is then subjected 
to a spatially uniform, time-dependent force. We are interested in the subsequent behavior 
of the system, particularly at strong driving. We avoid analysis in terms of quasienergies as 
is done in the Floquet method ||, as that method can only be used if the external forces 
are purely periodic. 

We solve the time-dependent Schrodinger equation numerically and monitor the prob- 
ability for the system to remain in its initial state as a function of time. We find that in 
regions of the parameter space characterized by strong driving, the system repeatedly ex- 
hibits a tendency to return to its initial state, leading to a series of recurrences in the initial 
state occupation probability; this effect is enhanced if the excitation frequency is above the 
threshold for ionization. These recurrences appear to have no classical analog. 

Quantum recurrences in the dynamics of wavepackets in time- independent situations 
have been extensively discussed in the chemical literature |§ , and have also been studied in 
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the context of quantum chaos 0. For time-dependent Hamiltonians, Hogg and Huberman 
H have established a recurrence theorem: a system described by a Hamiltonian with a 
discrete spectrum subject to a nonresonant, bounded and time-periodic potential will return 
arbitrarily close to its initial configuration infinitely often. In our case, the Hamiltonian has 
a partly discrete and a partly continuous spectrum and so the above theorem does not apply. 

We study a particle of mass m in a one dimensional potential under the influence of an 
external force, F(t). For technical reasons, we use a vector potential to represent the forcing 
term. The Hamiltonian is 



Ti = hujn 



2d 2 V 2 r . , , fx d 
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where we have introduced length and frequency scales x$ and loq with h 2 /2m = huoxl, 
and dropped an irrelevant time-dependent constant from the right hand side of Eq. ([[]). 
Hereafter, we work in dimensionless units and choose F(t) = Focos(ut). The Hamiltonian 
becomes 

where S = 2Fq/uj is a dimensionless measure of the strength of the forcing term and Q(t) = 
sm(u>t). 

The numerical method that we adopt to solve this problem begins with an expansion of 
the wavefunction in a set of time- dependent basis states rj k (t): 

*(t) = j2Mt)vk(t) (3) 

k 

In order to account exactly for the influence of the external force, we choose these basis 
states such that they evolve according to the field-dependent piece of the Hamiltonian 

<£-^(«)£. (4, 

The advantage of this procedure is that the choice of the basis at t = 0, ^(0), is dic- 
tated entirely by numerical convenience. We choose %(0) = exp(ikx), so that f]k{t) = 



77fc(0) exp(— iSjfyk) where j(t) = J" <?(£') (it'. These considerations lead to a set of equa- 
tions for the amplitudes X k {t): 



For our numerical studies, we choose a smoothly varying potential, V(x) = — Asech 2 (a;). 
For A = 0.75, this potential admits one bound state ipb(x) = tt^ 1 ^ 2 sech 1 / 2 (x), with energy 
Eb = —0.25. We integrate Eq. fl5|) numerically using the Crank-Nicholson algorithm || 
to ensure unitarity and monitor the probability to return to the ground state, P(t) 



For low strengths (S < 0.5) we find that our numerical results agree with the predictions 
of first-order perturbation theory. For example, for frequencies above threshold {u > 0.25), 
we find that the probability of remaining in the ground state decays exponentially with a 
decay rate given exactly by Fermi's Golden Rule. Closer examination of the decay curve 
shows that superimposed on the exponential decay is a small periodic ripple whose frequency 
is twice that of the external force. The amplitude of this modulation drops rapidly as the 
strength of the external force is decreased. This small periodic term arises from coupling 
between the continuum states induced by the external force and is dropped in the typical 
Golden Rule calculation when the "rotating wave" approximation is made. 

In the high strength regime the new phenomenon of recurrences appears. Figure [l], which 
shows the ground state occupation probability as a function of time for different values of the 
excitation frequency, clearly shows that perturbation theory fails to describe this situation. 
While we would intuitively expect rapid and complete decay, we see a surprising sequence of 
recurrences in the initial state occupation probability after nearly complete depletion. The 
recurrences become more pronounced at frequencies larger than the ionization threshold, 
although we expect that a particle supplied with sufficient energy to overcome the attractive 
potential should escape. Contrary to intuition, the system comes back to its initial state even 
though the density of sufficiently energetic quanta is substantial at these forcing strengths. 





where V(k) is the Fourier transform of V(x). 



(m\m)\ 2 - 
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We will now focus on the origin of these above threshold recurrences. At S = 10, the 
potential is clearly overwhelmed by the driving term. The system behaves largely like a 
quantum free particle subjected to the forcing field, the potential being responsible only for 
minor corrections to that behavior. We therefore proceed by first solving the problem of the 
free particle in the forcing field exactly and then adding in the potential through a Born 
approximation. 

The momentum space wavefunction i[>f(k, t) of a free particle in the external field can be 
found from Eq. (fj) with V = as 

i>f(k, t) = #(k, 0) exp[-i(k 2 t + Sj(t)k)], (6) 

where the initial wavefunction, given by the Fourier transform of ipb(x,0), is ip(k,0) = 
(l/7rV2)|r(l/4 + iA:/2)| 2 . 

A first Born approximation to the correct solution (valid in the case when the potential 
is small compared to the driving term) is a solution ^ of 

i^-[& + SG(t)k]1, B = V1> f , (7) 

where the approximation lies in including only ipf instead of on the right hand side. 
We solve this linear differential equation for ipB numerically and evaluate the approximate 
ground state occupation probability. Figure [|(a) compares the complete solution with the 
solution obtained by letting the initial wavepacket evolve under the influence of the forcing 
field only, while Fig. 0(b) shows the improvement gained with the Born correction. We see 
that most of the behavior of the exact solution is reproduced by the free particle result 
\{ip(0)\ipf(t))\ 2 . The bimodal structures present in the full problem are missed, but the 
overall size of the recurrences is well reproduced. So the existence of decaying recurrences is 
primarily due to the dynamics of wavepacket evolution in the external field and is relatively 
independent of the potential. 

While the results we have shown here are for a potential with only one bound state, we 
have also studied the potential V(x) = —6 sech 2 (x), which has two bound states at Ej, = —4 
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and Eb = — 1. For this problem we see qualitatively similar recurrences at S = 10, although 
with more structure, for excitation frequencies of u — 3, 4, and 5.61. Note that only the last 
frequency is above threshold, and the first one is the resonant frequency for the transition 
between the two bound states. 

To study the dependence of the recurrences on the form of the initial state, we examine 
the behavior of a free particle under the influence of the external force for an initial state 
given by a Gaussian wavepacket, i[)q(x) °c exp(— x 2 /4a 2 + ik x), of width a and momentum 
ko- The recurrence probability Pc(t) is given by 

(Sj(t) + 2k t/a) ' 



Pa{t) = i exp 
k + t 2 /a 4 



4 + t 2 /a 4 



(8) 



In Fig. |] we show the behavior of Eq. (Q) for ko = 1 and a = 1. These results are 
quite similar to those of Fig. and clearly show that recurrences are generic features of 
strongly driven quantum systems. We also see from Eq. ([8]) that recurrences will appear 
when the impulse given the particle in one cycle of the external field S^t) dominates the 
characteristic momentum of the particle. Furthermore, the decay of the recurrence height 
is controlled by the prefactor in Eq. (^) which arises from spreading of the wavepacket. In 
general, as we can see from Eq. ([]), if the argument of the exponential is stationary for k 
values where ipf(k,0) is large, then the wave function has a large overlap with the initial 
state and the particle shows recurrences. 

One further question that naturally arises is whether these recurrences have a classical 



analog. This is particularly interesting since Davis and Heller |10j have studied wavepacket 
evolution in several model potentials in time-independent situations and conclude that 
there is indeed very good correspondence between the time average of the quantum overlap 
|( , 0(O)|'?/'(t))| 2 and its classical analog. Furthermore, a classical mechanism involving the re- 
turn of an electron to the origin which then scatters out a second electron has been invoked 



IlJ in a model for strong-field double ionization of atoms |T2j . 

To compare the classical and quantum dynamics, we consider [ ]TTJ] , |r5]] a cloud of classical 
particles whose initial conditions are chosen so that their probability distribution mimics 



\ipo\ 2 as much as possible. Since we have chosen an energy eigenfunction as an initial state, 
we have several choices for the initial distribution of the classical cloud. We studied two 
different classical distributions. In one, we fixed the position probability density to match 
|V>(;e)| 2 and determined the momenta from the fixed energy, giving positive and negative 
momenta equal weight. In the second case, we generated initial conditions with position 
and momentum probability densities corresponding to Gaussians; the width of the initial 
position distribution was matched with that of the initial quantum state. In either case, we 
could represent the analog of the quantum overlap by using either the overlap of the time 
evolved position probability density with the initial one, / dx p(x,0)p(x,t), or an overlap 



integral QTOJ in phase space, Tr [p(0) p(t)). 

We find that in all cases there are no similarities between the classical and quantum 
results; the classical distributions fail to show recurrences for the same values of the pa- 
rameters of the forcing term. We thus conclude that the recurrences discussed here are 
purely quantum mechanical in nature, and that the classical model for the recurrence used 
in studying the double-ionization of atoms |TT| dramatically underestimates the amplitude 
for an electron to return to its original bound state, which may explain the discrepancy 



between the classical model and the experimental results BT2 



To examine whether the phenomenon of recurrences in strong driving fields exists in 
higher dimensional phase spaces, we have also considered an effective 2D problem: the 
hydrogen atom in its ground state subject to linearly polarized radiation. We observe that 
recurrences in the ground state occupation probability do exist at strong driving above 
the ionization threshold, even though they are much smaller. This work shall be reported 
elsewhere. 

We also expect recurrences even when the external force is not purely periodic. The 
crucial parameter that determines whether the external force is strong is the coefficient of the 
momentum operator in the last term in Eq. ([]]), / F(t)dt/(h/x ), which can be interpreted 
as the impulse delivered to the system in one cycle in units of the characteristic momentum 
of a particle in the potential. For above threshold excitation, when this number exceeds the 



strength of the potential, the system will behave like a free quantum wavepacket and show 
all resulting effects described here. 

We have therefore shown that whenever a quantum system is subject to a uniform force 
varying sinusoidally in time with a strength parameter that overwhelms the potential, above 
threshold excitation will lead to a set of decaying recurrences in the initial state occupation 
probability. While the exact nature of the potential and the shape of the initial packet deter- 
mine the detailed behavior and the long time decay in the magnitude of the recurrences, the 
existence of the recurrences is relatively insensitive to these features of the problem. While 
qualitatively similar decaying recurrences have been observed ||[7j in time-independent situ- 
ations, the ones we describe are more startling; the quantum system repopulates the ground 
state even when more than enough energy is available from the external field for escape from 
the well. 

This research was supported by the National Science Foundation under Grant No. DMR 
9020310. 
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FIGURES 

FIG. 1. Probability of initial state occupation as a function of reduced time. T = 2tt/uj is the 
period of the forcing function. Note that lo/loq = 0.25 corresponds to threshold excitation. 

FIG. 2. Agreement between approximate and exact calculations for above threshold excitation. 
T = 2tt/uj is the period of the forcing function. The exact calculation is shown in both figures by 
the solid line, while the dashed line shows (a) the free particle result |(V>(O)|'0/(^))| 2 from Eq. (||) 
and (b) the Born approximation \(ip(0)\tpB (t)}\ 2 from Eq. (|7|). 

FIG. 3. Recurrences in the time evolution of a Gaussian wavepacket subject to an external 
field, according to Eq. (^|) with ko = a = 1 for two different driving freqencies. 
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